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Abstract 

We present an ab initio study of the (001) interfaces between two insulating perovskites, the 
polar LaAlC>3 and the non-polar SrTiOa. We observe an insulating-to-metallic transition above 
a critical LaA103 thickness. We explain that the high conductivity observed at the TiC^/LaO 
interface, and the lack of similar conductivity at the SrO/A102 interface, are inherent in the 
atomic geometry of the system. A large interfacial hopping matrix element between cations causes 
the formation of a bound electron state at the Ti02/LaO interface. This novel mechanism for the 
formation of interfacial bound states suggests a robust means for tuning conductivities at various 
oxide heterointerfaces. 
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A wide variety of novel physical phenomena are observed at the hetero-ep it axial interfaces 
between perovskite oxides such as high conductivity |3j, |4|, external field doping 
magnetism [?| and supercoductivity sj], all of which are lacking in the bulk constituents. 
Understanding the rich physics of these systems, and their potential applications, has been 
the focus of intensive research in the past decade. 



Recently, the discovery of many novel properties at the (001) epitaxial inter 



LaAlOs/SrTiOs (LAO/STO) has boosted considerable research 
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Since both LAO and STO are perovskite oxides, two types of interfaces can form along 



the (001) direction: TiC^/LaO (n-type) and SrO/A102 (p-type) {4], E|. While experiments 
have found that the p-type interface is insulating, a high-mobility quasi two-dimensional 
electron gas (Q2-DEG) is formed at the n-type interface |4J when the LAO thickness ex- 
ceeds a critical separation p, [lo|. The origin of the Q2-DEG is the key to understanding 



the novelty of the LAO/STO inter 
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Broadly, there are two representative schools of thought. The first suggests an intrinsic 



mechanism due to the polar discontinuity at the interface 0, [9} . Since the nominal charges 
of the LaO and AIO2 planes in (001) LAO are +1 and -1 (STO has neutral SrO and Ti02 
layers), there is a macroscropic electric field through the LAO. The voltage across the LAO is 
approximately proportional to its thickness and diverges with increasing thickness. To avoid 
the divergence, half an electron(hole) per two-dimensional unit cell must transfer across the 
n-type(p-type) interface, respectively {4], 9]. Nevertheless, this symmetric description can 
not account for the drastic difference in electronic properties of the two interfaces. The 



second school of thought posits that extrinsic oxygen vacancies generated by the 
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laser deposition create the dominating carriers responsible for the conductivity [9j, 
However, both the existence of a lower limit in the carrier density of fully oxidized n-type 



nsity < 



interfaces [25|] and the insulating property of p-type interfaces |4J, |9J strongly suggest that a 
more fundamental intrinsic asymmetry exists in the electronic reconstruction 26|] between 
these two interfaces. In this Letter we reveal a fundamental difference in the electronic 
reconstruction between p-type and n-type interfaces and propose that the atomic geometry 
accounts for the difference. 

We begin by considering a schematic of the local energy gap in Fig. 1. We start with the 
simplest rap-type geometry in which several LAO layers are sandwiched between two thick 
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slabs of STO with an n-type interface on the left and a p-type interface on the right. In 
both STO slabs, the local energy gap equals that of bulk STO (3.2 eV). For the LAO film 
it has the bulk LAO value (5.6 eV) except that the macroscopic electric field through the 
polar LAO lifts up its local energy gap, as illustrated in Fig. la. The energy gap of the 
entire interface system is thus determined by the difference between the conduction band 
edge (CBE) of STO on the left and the valence band edge (VBE) of STO on the right. 
Therefore, the LAO film serves to diminish the gap. As the LAO film thickens, the voltage 
across the LAO accrues until the VBE of STO on the right exceeds the CBE of STO on the 
left (bottom panel of Fig. la). At this point, the CBE on the left become more energically 
favorable than the VBE on the right and electrons transfer across the LAO. These transferred 
electrons reside on the Ti <i-orbitals of the CBE on the left while holes are left behind on 
the O p-orbitals of the VBE on the right. This picture also applies to the n-type interface 
(Fig. lb) except that the holes stay in the O p-orbitals of the LAO surface valence bands. 
For the p-type interface (Fig. lc), the polarity of the LAO film is reversed so the local band 
edges bend down, and the electron transfer occurs when the CBE of LAO is lower than the 
VBE of STO. In all three cases, the condition for electron transfer is that the voltage across 
the LAO film exceeds the band gap of STO (np-type and n-type) or the band gap of LAO 
(p-type). The minimal LAO thickness satisfying this condition is the critical separation. 
To verify and move beyond schematics, we employ density functional theory within the ab 



initio supercell plane- wave approach [27j with the local density approximation (LDA) [28], 
and ultrasoft pseudopotentials 29 ] . All the super cells 33 1 in our simulation have the format 
(STO) Sl / (LAO)// (STO) S2 /vacuum, in which s\, S2 are the numbers of STO unit cell and I is 
that of LAO. Si, s 2 could be half integers due to an extra atomic layer (i.e. half a unit cell). 
The vacuum separating periodic copies is ~ 15 A in the (001) direction. For the n-type 
and p-type interfaces, s 2 is set to zero. For the np-type, the interface between (ST0) S1 
and (LAO); is n-type and that between (LAO)/ and (STO) S2 is p-type. Typical values of 
Si and S2 are 5 (for np-type) and 11 (for n-type and p-type); I is varied from 1 to 7. The 
interfaces are along (001) so the z axis is perpendicular to the interfaces. The lengths in x 
and y directions of the unit cell, subject to periodic boundary conditions, are fixed to the 
theoretical lattice constant of STO a = 3.85A (1.5% smaller than the experimental value). 

n 

All atoms except two STO layers facing the surface [30] are fully relaxed in the z direction 
until every force component is smaller than 26 meV/A. The charge densities in Fig. 2 have 
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been checked with convergence by reducing the force components to 13 meV/A or smaller. 

Our ab initio calculations support the above schematic view. For n-type interfaces, 
once the transfer occurs, electrons accumulate at the interface and inhabit electronic states 
composed of Ti <i-orbitals in the STO; for a (001) growth direction, the states have clear 
d xy character. For p-type interfaces, holes inhabit states of pure O p-character in the STO. 
Electron transfer occurs at 3, 4 and 5 LAO unit cells (u.c.) within DFT-LDA for np-type, 
n-type, and p-type interfaces, respectively. Table 1 lists the energy gap A versus the number 
of LAO u.c. as well as the corresponding potential change V per added LAO u.c. Based on 
the data, V extrapolates to 0.7 eV. 

This potential change is predicted accurately by DFT as it is a ground-state property 
depending only on the electron density. However, since DFT describes excitations poorly, 
the computed STO and LAO band gaps are underestimated (1.85 and 3.25 eV, respectively) 
relative to experiments. Therefore, we use the experimental band gaps of STO and LAO 
to estimate the actual critical separation. The difference between the experimental and 
calculated band gaps is taken into account by adding extra LAO layers, each of which 
reduces the energy gap by 0.7 eV. Thus we predict experimental critical separations of 
5, 6 and 8 LAO u.c. for np-type, n-type and p-type interfaces, respectively 3JJ. The 
critical separation for the np-type interface agrees very well with experiment 10J while our 
prediction for the n-type interface (6 u.c.) is larger than the experimental value (4 u.c.) jlj]. 
An important physical reason might be that we simulate perfectly flat, periodic, and defect- 
free LAO surfaces exposed to vacuum. In experiment, however, unavoidable imperfections 
of the surfaces (e.g. due to reconstructions or defects) may be responsible for an abrupt 
transition in the conductivities at smaller critical separation. For example, Cen et al. report 
that the creation and annihilation of oxygen vacancies at the surface could change the polar 
structure of the LAO film and the electronic properties of the interface Ijj]. 

Our results also provide clear evidence of an asymmetry in the electronic reconstruction at 
n-type and p-type interfaces. Fig. 2a) and b) illustrate the x|/-averaged spatial distribution 
of the transferred electron and hole densities along the (001) direction 32]. At the n-type 
interface, the electron density decays monotonically into the STO, indicating a bound state. 
However, the holes are not bound to the p-type interface because their density maximum 
is located in the middle of the STO substrate. (This feature persists as we thicken the 
STO.) This suggests that holes diffuse into the STO substrate while electrons stay close to 
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the interface. Because carrier mobility depends on trapping and scattering, the asymmetry 
in these distributions provides a tentative explanation for the observed difference in con- 
ductivity at these two interfaces. Under an applied field in the xy-plane, electrons scatter 
only from defects within a few nanometers of the interface but the holes scatter off defects 
throughout the entire substrate. Given the high quality of the interface [9], electrons should 
be much more mobile than the holes. Furthermore, holes will easily be trapped on inevitable 
oxygen vacancies in the substrate which will further reduce their mobility. 

There is a single experiment that infers the presence of 32 ± 6% oxygen vacancies per 
interfacial unit cell and attributes the insulating behavior of the p-type interface to the 
vacancies j^]. However, our DFT calculations show that oxygen vacancies, just like holes, 



prefer to stay away from the p-type interface 



33| . Thus simulating a clean p-type interface is 



certainly consistent theoretically. The discrepancy between experiment and theory needs to 
be resolved by additional experiments and theoretical calculations. Nevertheless, we agree 
that oxygen vacancies tend to trap and immobilize the hole carriers and thus render the 
p-type interface insulating. 

Then this begs the question: why are the electrons bound while the holes are not? 
Our answer hinges on topological asymmetry. The holes reside on O p-orbitals and the 
oxygen atoms form a continuous network across the LAO/STO interface (Fig. 2d). We 
have calculated the hopping matrix elements between neighboring O p-orbitals and found 
no significant discontinuity at the interface. Hence, the voltage across the polar LAO drives 
holes into the STO, but once there, they delocalize into the substrate. 

However, the situation is very different for electrons. In STO and LAO, the conduction 
band edge states reside on the Ti and La cations, respectively. The Ti occupies the B site 
in STO while in LAO the La occupies the A site (in the standard ABO3 nomenclature 
for perovskites), resulting in a discontinuity in the arrangement of the electronically active 
transition metals across the n-type interface (Fig. 2c). A detailed examination shows that 
the electrons occupy mainly the Ti rf xy -orbitals but have a small weight on the La d xy - 
orbitals in the first LaO layer at the interface. We calculate that the interfacial hopping 
matrix element t between Ti (i^-orbitals in the Ti02 plane and the La d^-orbitals in the 
LaO layer (Fig. 2c), which does not exist in either bulk oxide and is unique to the interface, 
is 0.7 eV, about 100 times larger than that between neighboring Ti d^-orbitals in STO. 
Two factors make this hopping so large: first, La 5<i-orbitals are spatially more extended 
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than Ti 3d-orbitals, and second, the distance between the La and Ti is 86%(\/3/2) of the 
Ti-Ti distance in bulk STO. 

Thus, as a zeroth order approximation, we focus only on the interacting Ti and La at the 
interface. The strong interfacial hopping causes the Ti d xy state to lower its energy and the 
La d xy energy to rise. In chemical language, we form a pair of bonding and anti-bonding 
states between the two cations across the interface. Therefore, compared to the other Ti 
(i-orbitals in the substrate, the Ti d-orbitals at the n-type interface are more energetically 
favorable and bind electrons to the interface. 

First principle simulations provide complete control over the system and allow us to 
verify the above picture via a gedanken experiment. We take the optimized (relaxed) n-type 
interface structure as the reference. We manually move only the Ti atom at the interface 
by an amount z, allow the electronic system to achieve self-consistency, and calculate the 
new interfacial properties including the hopping matrix element t. When z > 0, the Ti 
atom moves away from the interface, the Ti-La bond is elongated and t is weakened. When 
z < 0, the Ti atom moves towards the interface, shortens the La-Ti bond and strengthens 
t. Fig. 3 shows how the transferred electron density responds to changing z. As the Ti 
moves toward the interface, the electron density further localizes on the first Ti atom at the 
interface. When the Ti moves away, the biggest peak in the density profile then appears on 
the next Ti atom away from the interface. These results provide strong direct evidence that 
the large interfacial hopping between Ti d- and La d-orbitals traps and binds the electrons 
close to the interface. 

In conclusion, we reveal a fundamental asymmetry in the electronic reconstruction at 
the n-type and p-type LAO/STO interface, which is accounted for by a strongly enhanced 
electronic hopping between the interfacial cations that host the carriers. This mechanism 
opens a pathway to nanoengineering electronic states at heterointerfaces. For example, 
doping different cations at the n-type interface would change the spatial extent of atomic 
orbitals, modify the hopping element and thus tune the binding of the electrons to the 
interface. Mechanically, cation motions relative to the interface, e.g. by a ferroelectric or 
pieozoelectric effect, could similarly modify the distribution of the interfacial electron states. 

We are grateful to Charles H. Ahn, Frederick J. Walker, Victor E. Henrich, Eric I. Altman 
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FIG. 1: (Color online) Schematic local energy gaps of a) np-type, b) n-type and c) p-type interfaces, 
respectively. For b) and c), the righthand LAO surfaces face the vacuum. The shaded parts are 
filled valence bands and the upper parts are empty conduction bands. In a) from top to bottom, 
the energy gap of the interface system A (highlighted in dashed purple/dark gray) is diminished by 
adding LAO layers until it disappears. For p-type interfaces, the energy gap is not reduced until 
the conduction band edge of LAO becomes lower than that of STO. a is the STO lattice constant. 
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TABLE I: DFT-LDA computed energy gaps and potential changes versus the number of LAO layers 
i for np-type, n-type and p-type interfaces (i=0 refers to a pure STO substrate). A is the energy 
gap of the interface system with the corresponding number of LAO layers. V* is the macroscopic 
potential change due to adding the i-th LAO layer. 
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FIG. 2: (Color online) a) The total transferred electron density integrated over the xy-plane in the 
STO conduction bands at the n-type interface. The electrons are bound to the n-type interface and 
decay into the STO substrate, b) The total transferred hole density integrated over the xy-plane 
in the STO valence bands at the p-type interface. The holes are delocalized with the maximum 
around the center of STO slab. Interfaces are highlighted by the colored dashed lines. Layers are 
measured in units of the STO lattice constant a. The integrals of these charge densities over the 
STO region are normalized to unity, c) Schematic of different types of hopping between neighboring 
Ti ci-orbitals and between Ti d- and La (i-orbitals at the n-type interface, d) Schematic of hopping 
between different O p-orbitals at the p-type interface. 
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FIG. 3: Evolution of the total transferred xy-integrated electron density in the STO slab and the 
interfacial hopping matrix element i as a function of the movement of the Ti atom at the n-type 
interface, z > 0: Ti moves away from the interface; z < 0: Ti moves towards the interface. The 
integral of each electron density over the STO region is normalized to unity. (Distances and the 
movement z are measured in units of the STO lattice constant a.) 
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